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Abstract 



Qh Spin interactions in relativistic quark-antiquark system in mag- 



netic field is considered in the framework of the relativistic Hamil- 
tonian, derived from the QCD path integral. The formalism allows 
to separate spin-dependent terms from the basic spin-independent 
£> . interaction, contained in the Wilson loop, and producing confining 

"^ | and gluon exchange interaction. As a result one obtains relativistic 

~f~, ■ spin-spin interaction V ss , generalizing its nonrelativistic analog. It 

O ■ is shown, that in large magnetic field eB, V ss modifies and produces 

t^j- \ hyperfine shifts which grow linearly with eB and preclude the use of 

perturbation theory. We also show, that tensor forces for eB ^ are 
active in all meson states, but do not grow with eB 



1 Introduction 



The spin-dependent terms in the interaction of two fermions in QED have 
a long history [H [21 El H] , where spin terms were deduced based on Pauli or 
Dirac form of Hamiltonians. In particular the hyperfine (hf ) interaction term 
was introduced by Fermi in [3] in the nonrelativistic form. Later on the QED 
perturbation theory was used to consider hf interaction in higher orders j6] 
see [7] for a review, and also in presence of magnetic field (m.f.) [8]. 

An important tool for the higher order relativistic treatment of bound 
states in QED was Bethe-Salpeter equation (BSE) [9], or its instantaneous 
Salpeter form [TO] . 



It is the purpose of the present paper to consider spin- dependent effects in 
the dynamics of the quark-antiquark (qq) mesons in strong m.f. The physical 
interest of such systems is due to high magnetic fields, which are expected in 
neutron stars, early universe and heavy-ion collisions, see [11] for discussion 
and references. 

In such QCD systems, as mesons or baryons, the main part of dynam- 
ics is of nonperturbative (np) origin, e.g. the confinement interaction V con {, 
which cannot be described by perturbative propagator particle exchanges. 
In addition the relative time problem existing in BSE is more fundamental 
in QCD, where the vacuum field correlations with period A, yielding confine- 
ment (and all strong interactions, including spin-dependent ones) are shorter 
in time, than physical particle exchanges, so that each quark in meson prop- 
agates in the already time-averaged vacuum [12]. Therefore to establish the 
fine and hf structure of strong interactions, one must use instead of BSE 
another relativistic formalism, namely the relativistic path integrals method, 
based on the so-called Fock-Feynman-Schwinger Representation [131 IE]- Re- 
cently a new form of this method was derived in [T5] , where a new integral 
representation of the qq Green's function was given. 

Basically, the main part of the strong dynamics (non spin-dependent) is 
given by the Wilson loop average and the methods, based on FFSR, give 
all meson, baryon, glueball and hybrid properties in good agreement with 
experimental data, see [TB], [T7] for reviews. 

Of special importance for us here are spin-dependent forces in QCD 
(which we shall test applying high m.f.). The latter have been derived for 
large quark masses in [18], and in the framework of our method for any 
masses in [El [191 EH EE1 [221 123] . 

The most important development in the study of the spin-dependent (SD) 
interactions consists of three advancements: 

1. Definition of perturbative and np, SD contributions and expressions for 
them in terms of the standard field correlators. This was done for large 
masses in 



2. Derivation of SD terms in the relativistic systems for low and zero mass 
quarks, done in [TBI [2U] and 



3. Definition of SD terms for nonzero temperature and in the deconfined 
QCD phase, which was done in 



4. Reexamining of the SD forces and suppression of spin-orbit forces in 
hadrons by the string motion, [23J. 

In the present paper we perform calculation of SD forces in the relativistic 
qq system in the presence of m.f. eB, which strongly modifies these forces, 
and in particular yields a linear growth in eB for the hf term. We also show, 
that m.f. creates deformation of the meson shape and thus stipulates tensor 
forces in the originally spherical meson. 

A possible contradiction of these effects with the positivity property of 
the Green's function is formulated and possible outcomes are discussed. 

In the same formalism we calculate the selfenergy correction to the meson 
mass in strong m.f. 

The paper is organized as follows: in section 2 we define spectral proper- 
ties of Green's function in m.f., in section 3 we derive the qq Green's function 
and relativistic Hamiltonians in m.f. from FFSR, in section 4 we concentrate 
on SD forces in m.f. and write down explicit expressions for them. The sec- 
tion 5 is devoted to tensor forces in m.f. The last section contains discussion 
of the results and an outlook. One appendix contains details of derivations. 

2 Absence of pair creation in color Euclidean 
and magnetic fields 

Consider quark-antiquark system in the external magnetic field, interacting 
with color vacuum Euclidean and perturbative fields, so that the covariant 
derivative is 

D, = D^Af, A^t a ) = d, - ieAf{x) - igA^t", (1) 

where Ay can be decomposed into the magnetic ("Euclidean") part A^f' 
and the electric part, e.g. 

4 6) =4 B) + 4 e) (pert), 

A (B) = i(xxB), 4 B) = (2) 

for the constant magnetic field B along z axis. Then the partition function 
averaged over nonperturbative vacuum can be written as 

(Z) A = ( I D^D^P exp(-#(m + D)$)) A , (3) 



where the averaging over gluonic fields is implied, 

(K) A = ji^exp (-1 j(F;^)d 4 x) DA, (4) 

and we disregard for simplicity gauge fixing and ghost terms. 

Integration over D^pDij) inQ yields the standard answer (where the 
proper renormalization is implied), 

(Z) A = /exphrln(m 2 -D 2 )\ = /exp (Ur J°° — e ~<™ 2 -&)\ \ . (5) 

Now the question of stability of the vacuum in the given external fields 
can be associated with the nonnegativity of the operator (m 2 — D 2 ), since 
otherwise negative eigenvalues of this operator would provide imaginary part 
in the exponent of (jSJ), which implies finite probability of pair creation, as it is 
clearly seen in the Schwinger expression for the pair creation in the constant 
electric (nonEuclidean) field [24] . 

In what follows we show that in purely Euclidean (colorelectric or col- 
ormagnetic) fields and in the magnetic field (external and perturbative) the 
operator (m 2 — D 2 ) is nonnegative and hence pair creation is absent and 
vacuum stability is ensured. 

To this end consider the Euclidean operator 

iD^ = 7m (»0„ + eA$(x) + gAp a ), (6) 

with Euclidean and hermitian 7 matrices, 74 = /5; 7$ = —i{3ai] 7+ = 7$, and 
with hermitian A^ and A A = A%t a ,A£ = At, while A^ = 0, so that 

(tD^)+ = iD^^ (7) 

and for the eigenf unctions u n and eigenvalues A n we have 

iD^Un = X n u n , \ n real. (8) 

Hence 

-b 2 u n = (iD^ tl ) 2 u n = \ 2 n u n , A 2 >0, (9) 

and || m 2 — D 2 \\> m 2 , which implies vacuum stability and no pair creation 
for any m > 0. 



In terms of the quark propagator S = jr 1 — , one has representations 
. c _^ u n (x)u+{y) ( 1 \ _ ^ u n (x)u+{y) 




Writing u n = , one obtains the following equations 



-D 2 + a{gE -gU- eB){ip n - Xn ) = A r 2 t (</? n - Xn), 

-D 2 - <r(gE + gH + eB)(<p n + Xn) = \ 2 n {<Pn + Xn), (11) 

where D 2 = (d a — ieA^ — igA^) 2 , and -A^, E, H correspond to the color 
fields e.g. A^ = A^t^p, while Ay is the electromagnetic field, and B is the 
external magnetic field. From the system (ITT]) one can see, that for hermitian 
Euclidean fields A^\ A n , E, H, B the eigenvalues A 2 are real and the Green's 
function 2 ]_ ~ 2 has only positive set of eigenvalues. 

Note, that in real electric field A Q , Af = iA Q and the property of 
nonnegativity of (-D 2 ) is violated, implying possible quark pair creation, 
while in absence of Aq , but with real vacuum color field A\ = v4 4 and 
hence real colorelectric vacuum field EJ ac ~ d%A^ vacuum is stable with 
known vacuum condensate ((EJ ac ) 2 + (HJ ac ) 2 ) . However, the quasizero modes 
with small X n can accumulate, implying Chiral Symmetry Breaking (CSB), 
signalled by the Banks- Casher formula [25], see [26] for review and discussion 
of CSB from this point of view. As we shall see, perturbative colormagnetic 
interactions in the lowest order, in the external magnetic field can violate 
positivity condition (9), signalling the divergence of the perturbative series. 

3 Relativistic qq Green's function in magnetic 
field in the path-integral form 

The quark Green's function in magnetic field can be written as 

S q (x, y) = (m + D)- l y = (m - D) x (m 2 - D 2 )^, (12) 

where 

m 2 - D 2 = m 2 - D\ - ga^F av - ea^FJfJ (13) 



and D^ is given in §§§, while 



^F„„=(<^ »?A „ F.(f> ^ B ° 



'/i^- 1 - /zl/ 



(14) 



v erE <xH J ' ^ ^ ~ ^ crB J 

One can use for S q the path-integral form [131 US] 

/•oo 

5,(a;, y) = (m- D) x / d Sj D 4 z$f > (x, y)e~ K (15) 

JO 

where 

$i F) {x, y) = P A exp (ie / A^dzi + ig A^dz^ + ea^ u / F$dr + ga^ / F^dr 

V Jy Jy Jo JO 

(16) 



if 



2 . 1 /^V 

m+ 4U7 



dr. 



(17) 



Then the gg Green's function can be written as 
G qiq2 (x,y)=j o d Sl j o ds 2 (D 4 z^) xy (D 4 z^) xy (TW a (A)) A W (T (A^)) 

where 



(18) 



W a (A) = P exp(ig A^(z)dz^ + g cr^F^r^dn - g (T^ u F^(T 2 )dT 2 ), 

(19) 



W^M) = exp^ / Afdz^-ieJ A^dz^+eJ dr 1 {a, u F^)-e 2 dr 2 (a, u F$)) 
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T = -tr{Y 1 (m 1 - D 1 )T 2 (m 2 - D 2 )) exp(-X 1 - K 2 ) 



(20) 

(21) 



and T 1 = j^, T 2 = j u for vector currents, while 1^ = T 2 = 75 for pseu- 
doscalars, and the symbol tr implies summation over color and Dirac indices 
and refers to all terms. 

At this point we introduce new variables Ui in the path integral (TT5|) . 
defined via the connection between the proper time n and the real Euclidean 



time tf = z^iji) (see details in Appendix 1 of [15]), and integrating over time 
fluctuations using 

T _Tdui _ dtf 

2uj 2ojf: 2ujj 



roc roo At,]. 

Jo Jo 2uf 

(22) 
where ( ) a 24 means the averaging over time fluctuations, which can be written 
in terms of the averaged Wilson line 



and K{uj) is 



<<l>-.(.'--//)>A-., = \l^r**(x,y), (23) 






T = \x 4 -y 4 \. 

In this way the path integral in dsiD A z^> is replaced by du)i(D 3 z^'), and the 
qq Green's function ( TT8l) can be written as 

Z7T JO Wl ' JO W2 ' 

(25) 

where 

(TW) = 4irF(WW)"exp(-X(a;i) - Ar(w 2 )) (26) 

PU is the averaged over time fluctuations Wilson loop, see [15], and Y" = 
\Ti(mi — iPi)T2(m,2 — ipz). Note, that the operator ordering is not taken 
into account in (|26|) and we shall take care of it below computing SD terms. 
Using cluster expansion and nonabelian Stokes theorem [27] one can 
rewrite W„ as 



(TrW(C)) = (TrexpigJdTr^(z)F^(z)} 



(27) 
00 (ia) n r r 

^pE 1 ^ J d7r(l)...jdTr(n){(F(l)...F(n))), 



71=1 



where d-K pv = ds M „ + afydri — aj 2 Jdr 2 , and ds^ is an area element of the 
minimal surface, which can be constructed using straight lines, connecting 
the points zjp(t) and z^'(t) on the paths of q± and q 2 at the same time t 

[T2l [28] . Note, that z\ (t) = z\ (t) = t. Then the spin-independent part of 
the exponent reduces to the confinement term V con f(r) plus color Coulomb 
potential Vcoui, while spin-dependent part Vsd depends also on proper time 
variables ti,t 2 , (see (T2j [19] for derivation and discussion). For the case of 
zero quark orbital momenta with the minimal surface, discussed above, one 
obtains a simple answer for (W a (A)) a, which we shall derive below. 

The average ((. . .)) stands for connected correlators, for example, for the 
bilocal correlator, ((F(1)F(2))) = (F(1)F(2)) - (F(1))(F(2)), and F pv = 
d p A v — d u A p — ig[A p ,A u ] is the vacuum field strength. Obviously, due to 
the 0(4) rotational invariance and colour neutrality of the vacuum, ((F)) = 
(F)=0. 

In the Gaussian approximation for the vacuum, when only the lowest, 
bilocal correlator is retained, one has, with the accuracy of a few per cent 
(see Ref . [2H E3] for the discussion) : 



(TrW{C)) oc exp 
where 



-- / dn pu (x)d7c Xp (x')D puXp (x - x') 



(28) 



a 2 
D puXp (x - x') = j^((TrF^(x)^{x, x')F Xp (x')^(x' , x))). (29) 

This bilocal correlator of gluonic fields can be expressed through only two 
gauge-invariant scalar functions D{u) and Di{u) as [30] 



D puXp {u) = {^\8 vp -8 tJ , p 8 v x)D(u)+- 



d 

{u\8 vp - u p 5 X u) + 



du p 




DM 



(30) 

The correlator D(u) = D(uq, |u|) contains a nonperturbative part and it is 
responsible for confinement: the QCD string formation at large interquark 
separations. Indeed J ds i4i (u) j dsa{v)D{u—v) = S J d 2 (u—v)D(u—v), where 
5* is the total area between the averaged quark trajectories. The fundamental 
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string tension can be calculated from the area law asymptotics of (J2SJ) for 
large area loop and is expressed as a double integral: 



poo /"OO 

d 2 {u-v)D{u-v) = 2 du d\D{u,X). (31) 

JO Jo 



Using (fT9l) with cr^F^ = 0, one obtains spin-independent terms in the 
qq interaction [3D] . 

V (r) = V eonf (r) + V OG E(r), (32) 



pT pOO 

V conf (r) = 2r dX dvD(X, v) -> or, (r -> oo) (33) 

jo jo 

^oge = f ArfA r d^r 4 (A, i/) = ~-, (34) 

jo jo 3 r 



where we keep only perturbative part of D\ and to the lowest order D\ er (A, v) 
We now turn to the spin-dependent terms of the qq interaction, 



16a s 



3tt(A 2 +^ 2 ) 2 ' 

Vsdi an d we shall be interested only in the zero orbital moment states for 
simplicity, hence only spin-spin interaction term V ss will be treated below. 

The spin-dependent terms in the interquark interaction are generated by 
the combination a^ v F^ v present in Eg. (fill) and therefore one needs correlators 
of the colour-electric and colour-magnetic fields, as well as mixed terms, 
separately They immediately follow from the general expression ( 1301) and 
read [31] : 

9 ■({TrE i (x)QE j (y)&)) = <% D E (u) + D E (u) + u\— \ + u iUj ' 



2 



n 2 ( r)D H \ 8D H 



^((Trtf^i^c&t)) = Sij \D H (u) + D?(u) + u 2 ^J - ^ ^ ;! 

^((Tri^x)^^}} = eiifcW4Ufc -^L_, (35) 

where u^ = x^ — y^, u 2 = u^u^. We keep here the superscripts E and if in 
the correlators D and D% in order to distinguish in principle the electric and 
magnetic parts of the correlators and thus to be able to consider a nonzero 
temperature T and to distinguish Euclidean and Minkowskian contributions. 
Indeed, while D E = D H and D E = D H at T = 0, at higher temperatures 
they behave differently. In particular, above the deconfinement temperature, 
T > T c , the electric correlator D E disappears, whereas Df and the magnetic 
correlators survive. 



It is clear from (Q3J), that in the norelativistic limit, when in (II 7p in 
/ cr^F^dr, dr { = ^-, only the upper left corner of (TH1) . i.e. (crH), will 
contribute to V ss . The corresponding derivation was done in [T2| [19] and 
gives 

T/ , s cr x tr 2 3(o-ir)(o- 2 r) - <Ti<T 2 r 2 

M r = 775 w) + — = F 3 (r) 36 

12mim 2 \lm\va2r 1 

where 

V ^ ) = -L dur dr> "' (37) 



*(r) = jT A, (3DT(r, ,) + ^^^ ) = ^^(r), (38) 

It is our purpose below to calculate V ss in the relativistic qq system and 
in arbitrarily large magnetic field B, and to this end we shall use below first 
the path integral formalism [TU [15], deriving the general structure of the qq 
Green's function, and in the next section we shall determine how relativistic 
V ss expressions depend on B. 

First we need to find the the Euclidean action S E g 2 in terms of 001,002 and 
common time t E of the qq system at if = if = t E . To this end we define the 

Euclidean Lagrangian L*^. We write =j—- = 2ui~^- = 2u>iZk, k = 1,2,3. 
Then all terms in the exponents in (TT9I) . ( 120]) and (I2TJ) can be represented as 
exp(— / dt E L E li}2 ) and thus we arrive at the following action 



/•Te 

S qi& = J dt 



^ + E(f(^) 2 



+sl 



ie t A k z k +2 ^ + 2 ^+ ei ^ +e 2 2 ^ +a|z z | g ^ _ ^ 

- (39) 

where S E contains (crF) terms. Here A k is the fc-th component of the 
QED vector potential, a is the QCD string tension and the contribution of 
terms (cr\F), (o" 2 -F) is separated in S E . The next step is the transition to 
the Minkowski metric and the construction of the Hamiltonian. This is easy, 
since confinement is already expressed in terms of the string tension. We 
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have exp(— / L E dt E ) —¥ exp(i J L M dtM), ^e —* Hm, an d 

dL M 
I 
As a result one obtains (back in the Euclidean time T — \x± — y±\), 
T f°° dwi r°° doj 2 

3/2 / ~ 

H M2 = H + H a + W, (42) 

W = V con{ + Vqge + AM SE + AM SS , (43) 



H qiq2 = E ^vf - L M , P f = j~ W = <*$ + e^ . (40) 

i OZh. 



G ^ = hi:^i:^YWe-^\ y) , (41) 



where 



oj? 



2 ( p B-|(Bxz<'))) +m?+u, 
#o = E^ o-^ l L , (44) 

i=l ZUJ i 

2wi 2w 2 v ; 

Here the terms AM$e and AM SS are produced by Sj?, and we shall find 
them as a first order correction. But before that we must treat the w» de- 
pendence either in the path integral (125]) . or the Hamiltonian (|4"2"|) . In the 
path integral Wj play the role of quark energy parameters, and one can use 
the spectral decomposition in (|4"T|) to rewrite it as 



<?(*,!/) = f f^ r-^f;4*ry(x|n>e-^<n|y) J (46) 

At large T one can use the stationary point method, and one defines up- 
front the extremum values of M n (ui,u)2), namely for the Hamiltonian H, 

H + H a + V conf + V OGE = H; H* = mV>*, (47) 

and u\ is defined from the condition 

W^>| „=„, , = 1,2. (48) 
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To have an idea of the possible meson masses and the values of u\ , which 



we shall use below, it is instructive to consider as in [12] the main part of 
the Hamiltonian, i.e. 

H = H + H (T + V coni , H4, = M$ (49) 

and replace V con { by the quadratic term, 

KoBf(r) = or ->■ V conl (r) = | I — + 7 J . (50) 

where 7 is the variational parameter, yielding some 5% accuracy in the re- 
placement ( l50l . Then the resulting mass M can be found explicitly as 



M(wi, W2, 7) = ^.n, + ^ + ^ , (51) 



where 



1 

e n± ,n z - — 



/ e 2£2 + (2n ± + 1) + J [n z + - 



70" 

+ T> (52) 



As a result one can estimate the masses M and u)\ at large m.f., since the 
basic pattern is defined by relative signs of eB terms in £ n± ,n z and H ff . 
Indeed, for eB > a one can write M = £ i=1 2 ^+^+1^1-^3 and 



2w, 



Co' 



(0) 



fi i = */m? + |e i B|-e i <r i B + 0(- N /ff), M = ^i + fl 2 - (53) 



Thus, for the neutral meson with e 2 = — ex, and o"i 2 ,cr 2z = (++),u; ++ 
is growing as \fe~\B, while for the (H — ) state cu + _ is tending to a constant. 
One can also find the wave function 



fa) = . ; = exp ( -^ - ^ ] , r, = z« - z( 2 ), (54) 



1 / vl v 2 

^=^= exp ^ 

7r 3/2 r 2 ro y 2r ± 2r5 

and 

At large eB ^> cr, one has r\ m -^, r$ m const ~ -4=, and hence 



l^°)| 2 ^^7^> (eBXr). (56) 
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This is the focussing effect of m.f., which is most important in SD forces, 
as well as in other processes [33] . 

It is important to stress, that we have kept in H only those terms, which 
are the main part of interaction, and therefore in M^ and a/ -* they are 
treated to all orders, i.e. exactly. However, the terms V ss and AMse are 
considered only as a perturbation, and therefore one should substitute there 
the values u\ obtained from ( HTj) . ( HSj) . where V ss , AM S e do not enter. The 
Hamiltonians H are considered in [15], and below we shall derive both V ss 
and AMse in the relativistic qq system. 

4 The quark- ant iquark spin-dependent inter- 
action in strong magnetic field 

The advantage of representation f fl8|) . ([25]) lies in the fact, that the only 
place, where the Dirac 7 matrices enter, is the local term (m — D), and it 
can be assembled in the factor Y with due care, while all the rest nontrivial 
spin dependence is contained in the (crF) and {crB) factors (TH|) . We shall 
demonstrate below in this section, that the correlators (a^'F)(cr^ k 'F) with 
i = k define AMse, while those with i 7^ k, define V ss . 

Consider the Taylor expansion in powers of the color spin interaction 
gcr^F^ = g(aF), m 2 - D 2 = m 2 - D 2 - g(aF), D 2 = D 2 - e(^), 



1 1 1 1^1 



2 



m 2 — D 2 m 2 — D 2 — goF m 2 — D 2 m 2 — D 2 m 2 — D 

1 ^9(vF)—^-g(aF)—^- + ... (57) 



m 2 — D 2 m 2 — D 2 m 2 — D 2 

One can define in (1571) the selfenergy correction to the mass, 

Am 2 (x,y) = -g(aF)J * ) g(aF) y . (58) 

\m 2 -D 2 J xy 

Am 2 = / d i (x-y)Am 2 {x,y)- 1 Am 2 (x,y) = -g 2 a i a k ((H i H k } + (E i E k }) xy Go(x,y) 

(59) 
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where 



while 



^(H^Hkiy) + Ei{x)E k {y)) = 2D{x - y)5 ik . (61) 

It was shown in [32] that in absence of magnetic field and in the limit of 
small m and small vacuum correlation length A, one can replace f 2^02 ) 
by the free propagator (47r)2( 1 a; _ y )2 , yielding 



XL) 



Am 2 = - — , (62) 

Ti- 



ll (ja»j) yieias tor tne total mass M£' , i\M n = 2^ 
zero mass q and g 



then the correction ( |58l) yields for the total mass M^ ^ , AM n = Y,i ~#, f° r 



Mf(u; )->Mf(u;o)-— . (63) 

7rw 

Note, that u; plays the role of the integration variables in (125]) . (HT|) and is 
defined from the stationary point condition (|48|) . where M^ ) does not include 
Am 2 . However for large m and small |aj — y\ < A one should multiply ( |62l 
with the coefficient rj{m\) < 1 calculated in |32j . 

Consider now the case of constant magnetic field B along z axis. One 
can calculate the effect of magnetic field on the selfenergy correction, to this 
end expand 

( 1 m f ^ m 2 -Dl + ea 3 B aF(G + + G.)aF 
{aF) m 2 -Dl-ea 3 B ^ F) = ^ (m 2 - D^) 2 - (eB) 2 ^ ~> 2 ' 

where G + /_ = (m 2 — D 2 ± eB) X y, and m 2 — _D 2 ps m 2 + |ej£?|. As was shown 
in Eq. (l53|) for large e-B 3> a one has a large effective mass oj*,oj 2 ~ m 2 + 
2\eiB\ in one of the Green's functions G + , G- and the corresponding Green's 
function will contribute -^-r)(eB), where r)(eB) = r](J2\eiB\ + m 2 X) is the 
coefficient, introduced in [32], e.g. 7/(0) = 1,7/(5 GeV) 2 ) = 0.03, see appendix 
of [32J for explicit expression. The final expression for AM n can therefore be 
written as 

Q 

AM n = --^-(l + V (eB)) } (65) 
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where Uq corresponds to the Green's function G of a given quark, when spin 
terms (crF) are absent. One can see, that at large eB the selfenergy correction 
numerator decreases approximately twice, as compared to eB — 0, m — 0, 
since r\[eB S> a) — > 0. 

We now turn to the (qq) Green's function and write it in the form 

G qi g 2 (x,y)d 3 (x -y) = tr([ri(mi - ip\)e' HT T 2 (m 2 - ijh)]), (66) 

where we have taken into account, that (d—igA) acting in (166]) on the Wilson 
loop, can be replaced by the momentum operator [55] . 

For the case, when H does not contain 7^ matrices, noncommuting with 
(rrii — ipi), one can rewrite ( 166]) as in [33], but now taking into account spin 
and isospin nonconservation in m.f., one must keep the possible eigenvalue 
dependence on the spin projections, 

/ G qiq2 (x, y)d\x -y) = 5> r ® e r ) u [ n Jr / (67) 

J n,v ^-'"n 

with £ 75 = £1 = 1, Ey = eR and , and the index 1/ denotes a specific 
polarization and charge component of quark and antiquark, e.g. v = 1 for 

(ttw), < h I component. As a result the quark decay constant of the 7^ 

state is 

r , M2 iv c (y r ,)|^)(o)| 2 

(JrJ = — : : H — , (68) 

where £„ occurs due to w$ integrations in (jSJ), see appendix of [H], and 

(Y ru ) = i*r(r v (mi - ipOHma - »&)). (69) 

In the case, when if contain spin-dependent terms, and in addition de- 
pends on magnetic field B, one should be more careful with the ordering of 
operators (crF), crB in H and projectors (mi — ipi), (m 2 — ip 2 ). 

Correspondingly, (m — D)cem be rewritten as 

{m q - D) x ->(mi-ipi), (m q --D) x = m 2 -ip 2 , (70) 

where p% = (iooi,p), p 2 = — (iu 2 , — p) and we take into account, that D^ix) 
acting on <& z (x, y) yields d^ — > ip^. 

. . ( mi + ui -icrp \ . , 

mi - ipi = m 1 + Wi7 4 - ?p7 = . (71) 

V zcr P mi — 0J\ J 
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At this point we have at least two possibilities for the relative ordering of 



.2 



factors (to? — D ) and (to § — D) in G qq -. We shall define this ordering as 
Gqq(x,y) = (tr[T l S q (x,y)T 2 Sg(y,x)]) A = 

= (MFiK - D) x (m\ - D 2 ):ir 2 (m q - D) y {m\ - D)£]) A (72) 

We shall show below, that this ordering yields correct results for spinless 
and spin-dependent parts, which can be checked in the nonrelativistic limit, 
whereas other orderings lead to wrong answers. Now 

( m 2 -u) 2 -icrp \ , . 

to 2 - zp 2 = to 2 - c 274 - z P7 = (^ ^ m2 + U2 )- (73) 

If is clear, that in the nonrelativistic situation, p <C to, u = y/p 2 + to 2 — >• 
to, the product (r(mi — ipi)r(m.2 — W2)) tends to the nonrelativistic projector 

(o™ 1 '" 2 ' o)' < 74 > 

We now turn to the spin-dependent terms and make expansion of to 2 — 
D 2 = to 2 — -D 2 — ecrS — g(crF) in powers of g(aF): 

(m 2 1 -D 2 )- 1 = A B + A B g(aF)A B , A B = (to 2 - D\ - e^B)- 1 (75) 



(to 2 - D 2 )- 1 = A B - A B g(aF)A B . i A B = (to 2 - D 2 t - e 2 aB)- 1 (76) 
and (j72p can be rewritten keeping only one spin-spin term as 

G q g(x,y) = G < q ^(x,y)-(tr[T(m 1 -ip 1 )A B g(a 1 F)A B T(m 2 -ip 2 )A B g((7 2 F)A B }} A . 

(77) 
one readily obtains for T = 75 

M = tr[y 5 (rn 1 -ip 1 ){a 1 F)^ 5 (m 2 -ip 2 ){cr 2 F)] = tr[(m 1 -ip 1 )(a 1 F)( y (m 2 -ip 2 )(a 2 F)) T ], 

(78) 

(m -it \(nF\ ( ( mi + u i) aH i ~ iopo'Ei. ( m i + Wi)aEi - zcrpo-H! \ 
11 % P l >^V - y ^crpcrHi + (toi - wi)<rEi, ^po-Ex + (mi - Wi)<tHi J ' 

(79) 
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(( -■" )( F )) T — ( * cr P cr ^ 2 + ( m2 + ^2)crH 2 , zcrpcrH 2 + (m 2 + w 2 )crE 2 

I (m 2 — o> 2 )crE 2 — icrpcrH 2 , (m 2 — o> 2 )crH 2 — zcrpcrE 2 

(80) 
Combining ( 1T9|) and ( 18"U|) . one obtains (Hj = H(xj), z = 1,2). 

M = tr a {crHicrH 2 (2mim 2 + 2wiu; 2 ) - 2(crHi)crp(crH 2 )crp) + 



+crEicrE 2 (2mim 2 — 2wiw 2 ) + 2(crEi)crp(crE 2 )crp) + 

+ crEi<TH 2 (-2z<rp(a;i + w 2 )) + crHi<rE 2 (2i<xp(u;i + u 2 ))} (81) 

In what follows we disregard first the terms, containing (<xp) or (erp)(erp). 
In the nonrelativistic limit, |p| — > 0,Ui — > rrii, one has 

(M) A = 8m 1 m 2 (H i {x 1 )H i (x 2 )) A . (82) 

Field correlators are expressed via two scalar correlators D(x\ — x 2 ) and 
Di(xi — x 2 ) as in (l5o^) . 

Comparing with the standard definition for the nonrelativistic hyperflne 
(hf ) term, one has 

_(l) rr (2) ^oo 2 

V hf = - v} H \r), v} H \r) = / efa/^r<fli(x)fli(y)>, ^ee x 4 - y 4 , 

12mim 2 J-oo N c 

(83) 
where 

u fl = x^,-y fJi , /i = l,2,3,4; r=|x-y| = |u|, z/ = w 4 . (84) 

One can separate perturbative part in Di 

D!(x) = Dr\x) + 0?(x), DT rt (x) = ^ + O(o2), (85) 

and define the potentials 

V} H \r)= T d^(H l (x)H l (y)) = 
= |_°° <^ Ud(t, v) + 3Di(r, i/) + 2r 2 ^ r 2 ' ^ = V 4 (D) (r)+V 4 (1) (r), (86) 
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v} E \r) = |_~ du (3D(r, v) + SD 1 (r, v) + (3v 2 + r 2 )^^-) = v[ D) -V L 

(87) 
As it is known [21J, the nonperturbative part of D\ and D yield much 
lower input in V} ' , the leading part is due to Df ert , i.e. V} . Inserting 
this, one obtains 

^^(^^(r), (88) 



4 (1) W. 



and the result (1831) . ( 1881) coincides with the known nonrelativistic limit. 
Now we turn to the relativistic case, Ui 3> rrii. First of all we note, that 

y(E) s ^) P ert (r) = Z^^(3) (r) . (gg) 

o 
Looking at (18T1) . one can see that in the relativistic case, when Ui ^> rrii, 
there is a cancellation in the spin-spin interaction, in the combination 

2(m 1 m 2 + Ut^VJ rt + 2(mim 2 - Wiu^)^ * = 

= 2W1W2 • 2V 4 (1) + 2mim 2 2y 4 (D) , (90) 

and multiplying this result with As, A# as in ( 1771) . in the B = case we 
obtain for A B ps ^t, A B = ^y, 



i 



^SfM^+s™}- (91) 

Here u)j = — L , and fij is defined in ( 1331) one can derive that Ui > u\ , 

e.g. in the nonrelativistic limit w,- -^m,, f^ — > m^ and also a), — > rrii. The 
same happens in relativistic case, when eB 3> a. In what follows for u>i 7^ o> 2 
it is implied, e.g., that d>+_ = u>i(H — )u>2(H — )• 

We now turn to the case of nonzero f? and now take into account the 
noncommutative (2 x 2) terms in H a and in V/,/, which we write in the total 
mass as 

M = M - /xiO"i z + ^2^22 + ao-io-2, (92) 

where /x i = g,a = ^(Vf). 

For 7r°,p°(s 2 = 0) states one obtains a standard mixing of (H — | and 
( — h I states of (<7i 2 , <t 2 ^|, with Vhf, where now for B 7^ we distinguish 

a n = (+- \acr l( T 2 \ + -) = —l-^W), (93) 

lzw_| 
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a 22 = (- + lateral - +) = T^""^), (94) 

2a 12 = 2a 21 = (- + |acr lC r 2 | + -) = — — (Vf }, (95) 

We also define Mn , M 22 as follows 
Mn = (M - G«i + A*z) - an)) w+ _ ; M 22 = (M + {fi x + ^ 2 ) - a 22 )) w _ + . (96) 
Finally, from det(M — E) = , we obtain the eigenvalues of M, 



Ex, 2 = \{M n + M 22 ) ± y( M22 2 Mu ) + 4a 12 a 21 . (97) 

Here w + _ and a>_ + correspond to the diagonal states of A B A B , i.e. 
for neutral mesons to the stationary values of oj^ in the states with spin 
projections (H — | and ( — h | respectively. It is clear, that at large eB, 
the values of u behave differently, i.e. o> + _ ~ const, while w_ + ~ \/aB, 
and M 22 ~ ve5, hence at large e.E> the nondiagonal part of M in ( 1971) is 
decreasing, and M tends to its diagonal eigenvalues. 

E x {eB -»- oo) -> Mu, £ 2 (e£ -> oo) ->■ M 22 . (98) 

One can notice, that for the (H — | state Mn contains at large eB the fast 
decreasing part, —an ~ — ^ 2 (0), and the latter is large in modulus, ip 2 (0) ~~ 
eB, hence leading to the negative mass E\. 

This happens already for only colormagnetic contribution (HiHk) to V} , 
while our conclusion in section 2 was, that eigenvalues of (mj — uf) are posi- 
tive, together with the total mass eigenvalues of the operator ( ^ A9 A, ). 

° or- ^m^ — Dfm^—D^' 

Hence we conclude, that the perturbation theory in V^f breaks down at large 
eB and one has to replace it with some modified form. However, as was 
understood already in ^\ [3] even at B = the perturbation theory with the 
potential V}g'(r) ~ c5^(r) is diverging since Vj^ for any c < ensures infi- 
nite number of bound states, which are physically irrelevant. Therefore one 
should in any case take into account the relativistic smearing of the hf inter- 
action, which appears due to the time integration in (138]) . which is taken in 
(13"8"j) along the straight line, instead of the complicated relativistic trajectory 
of the quark with time fluctuations, see [To]. The resulting smearing length 
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A > X con f, where A con f is the scale of D(x)., connected to the gluelump mass 
[55] . A con f ~ 0.1 -r 0.15 fm. On the lattice A > a, a is the lattice unit, 
a « 0.1 -hO.24 fm. therefore we replace V^'(r) by a smeared out version, e.g. 

P)(r) = (-j=\ e~^ 2 - V ss (r) = c~5^(r), M =is(H2) GeV. (99) 

Using the wave function ipiv) from (|54|) one obtains for (V ss ), in the PS 

meson 

, T % v cu 3 8%a s /-,„„x 

(K a > = ,c= -. 100 

For /i — )• oo one regains the original answer, (Vg S ) — >■ ct/> 2 (0). 

Since r\ = -^ (l + (feW) , r ±(e-B — >■ oo) — >■ 0, and (V^ s ) tends to a 
constant limit at large eB, preventing in this way breakdown of the vacuum 
due to vanishing of the meson mass. 

5 The tensor forces in magnetic field 

As was established in the previous section, the hf interaction, has the form 
(19T1) and in the m.f. the coefficients =-^- transform into - — = — where (ik), (i'k 1 ) 
(cri z ,(T2z) = (H — ) or ( — h), for S z = see Eqs. (j55]) - (l^]l . We now turn to 
the tensor forces and discuss, how they are transformed in the m.f. As was 
found in [161 EH] the total spin-dependent forces without m.f. can be written 
as 

(1) (2) 1 

Vf + VT = — V 4 (r) + - (3(ff«n)(^n)-^)«rP))K,(r) (101) 

IZ0J1OJ2 1Z0J1UJ2 

It is clear from f llOip . that for the spherically symmetric qq states the term 
Vz{r) is irrelevant, and therefore requires a nonzero angular momentum. 

The situation is drastically changing in the m.f., since in this case the 
form of the wave function is distorted as in the elongated ellipsoid. 

This fact implies the appearance of tensor forces in m.f. even in the 
ground qq state with the zero angular momentum. This situation was studied 
for the hydrogen atom case in m.f. in [34J in the nonrelativistic treatment of 
the tensor forces in the qq) system. 
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We are using again the correlator technic and the perturbative correlator 
D\ cr , in terms of which the tensor interaction can be written as in (1371) . 

The expectation value of the tensor term in the ground state with the 
wave function ip(p, z) fl54|) can be written as 



{V s S) = s , k / dupdpdzM3n . nk _ 8 ik )^ 2 {p, z){-r 2 ' 



QixuiiUJ2 J-oo r ' " ' ' dr 2 (r 2 + v 2 ) 2 

(102) 
Using the relations 



i i d 2 



9 i 9 i ' 9 i 9 i / W^JC 

n r z + v z + a r z + v l + a Jo 



(r 2 + z/ 2 ) 3 2 da 2 
and the qq wave function in the zeroth approximation 



•• a.„a_ 



(103) 



^ 2 (p,z) = 7Ve^ % JV = .^ - (104) 



7r 3/2 r 2 r 2 ' 



one can calculate the sum 



(105) 
Here if is 

9 f°° 2u 4 du 
K = r\ - 2 . 106 

Jo (u 2 + l) 2 (u 2 + r -^f/ 2 

Note, that we have omitted here for simplicity in (V4) ss ~ o"icr 2 , the 
additional terms, which appear in (l9Tj) . The integral K in f II 6 j) can be done 
explicitly, using the relation 

dU l :lnf P + v gZ^ V (107) 



o (w 2 +p 2 )(n 2 + g 2 ) 1 / 2 2p^p 2 — g 2 \P — Vp 2 ~ Q 2 J 

It is important, that X < r^, therefore since r^ — )■ for eS — > oo, see eq. 
f )55|) . the relative role of the tensor term in H105[) is diminishing with growing 
eB, and the absolute value of (V^ s ) depends only on the values of u)%,u)2- To 
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find these values, we again as in (1921) write the total mass, but now with the 
addition of the tensor term and for any quark charges 

M = M + a ( T l( T 2 + c^a lz + c i S + K-^ ) ^ 2) -^a^ + 2a^a^ ) ee M + h, 

(108) 
where 

8nas N, c« = -^, c^ = - 6 -^-, (109) 



9uiU2 2u\ 2uj- 

■ ta -./(i-il (no) 



QU\U2 4 

Moreover, one must distinguish the values of coefficients a, o, q in different 
spin projection states (<Ji z a 2z ) namely the values of o>j in a, b, Ci in different 
(crW, af>) states. E.g. for S z = +1(— 1) one must write UiU 2 — > tD+ + (^ 2 _), 
while for S z = 0, one has the matrix elements as in (I93I) - (I95I) . In a similar 
way for our general case with tensor forces in (11051) . one can write 

(+-\h\ + -)= -a u - 2&n + c { 3 - eg 5 , 

(- + \h\-+}= -a 22 -2b 22 + c%-c ( 2 2 l (HI) 

(+ ._ \h\ - +)= (- + \h\ + -)= 2(oi 2 - 612) = 2(o 2 i - 621), 

and all diagonal terms have the corresponding a>», e.g. <jj\ = u 2 = u± for 
an, On, C11 , and in nondiagonal matrix elements Q\Q 2 = uj + -U- + . 

From det(/i — E) — with elements in fillip one obtains two eigenvalues 
of the spin-dependent part Ei, E 2 . 

The resulting expressions for Ei, E 2 coincide with those in fl97|) . when one 
makes the following replacements 

an ->■ an + 2on, a 22 ->■ a 22 + 26 2 2, a i2 -> a i2 - b 12 , (112) 

and //1 — > —c- l >. 

Eqs. ( )5T|) . (|52|) contain a prescription for the values of cjy, z, j = +, — , 
entering in Mn, M22, which is valid also in the case of nonzero tensor forces. 



6 Conclusions and prospectives 

We have obtained explicit relativistic expressions for the spin-spin interaction 
terms in the qq system in the arbitrarily strong m.f. As a by-product we also 
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obtained in the same formalism expressions for the np self-energy corrections 
in m.f. These formulas are generalizations of the previously found expres- 
sions for the qq mesons in absence of m.f., see e.g. [16], and we also found 
corrections to those expressions, see Eq. (19T1) . where the terms 3^3- 14 (r) 
and ^V^ (r) are new. 

It is remarkable, how m.f. changes the spin-spin forces. First of all, 
the matrix element of the hf term 14 (r) ~ ^-'(r), in the strong m.f. is 
proportional to the i/j 2 (0) - the probability of coming together of q and q, 
which in strong m.f. grows as eB 

This effect is known in nonrelativistic case, where the hf term is (14/) ~ 

mm ~ a/m anc ^ was discussed recently in [33] for the hydrogen case. 
However, in this case Mi = M p is very large and (Vhf(hydr)) is small and 



r^j 



m 



2 



this growth of (Vhf) was thought to cease in relativistic limit, eB 

However, as we have shown here in the paper, the growth of (Vhf) in the 
relativistic case for eB — > 00 is possible whenever u 2 in the denominator of 
{Vhf) ~ ^r- does not grow with eB, which occurs for the qq states, where 
magnetic moment terms compensate the growth of the rest part of mass, 
e.g. for the < H — | states of the neutral mesons, like 7r°. Then hf terms 
yield negative (—3(14/)) contribution to the mass, linearly growing with eB 
in modulus, thus giving the absurd negative mass result. 

To disprove this result, we have given in Section 2 the proof, that the 
term 14 , which causes this problem, cannot generate negative mass in 
any m.f., and we are coming to the conclusion, that this discrepancy occurs 
due to the use of the perturbation theory for the hf term, proportional to 
the delta-function. Therefore we have suggested in section 4 the smearing 
procedure of the hf term, which takes into account the relativistic Zitterbe- 
wegung of quarks and should be supplemented with a rigorous procedure of 
the summation or estimation of the whole perturbative series. We stress, 
that this situation occurs solely due to m.f., which cannot, as shown in sec- 
tion 2, provide the pair creation and vacuum reconstruction, in contrast to 
the real Minkowskian electric fields, which are capable of the pair creation 
(as in the Schwinger famous formula [23]) and of the vacuum reconstruction 
with emission of positrons (as in superheavy atoms with Z > Z CTit ). 

Another interesting new effect is the appearance of the tensor force effects 
at nonzero m.f., which can be tested both in atoms and mesons. However, 
as we have shown above, at large m.f. the tensor force contribution does not 
grow with eB, unlike hf term, and is always smaller than the latter. 
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In this way we have completed the main part of the strong dynamics of 
mesons in magnetic field for the zero angular momentum, started in [To] . 

The author is grateful for useful discussions to M.A.Andreichikov, S.I.Godunov, 
B.O.Kerbikov, V.D.Orlovsky, and M.I.Vysotsky and to the members of the- 
oretical seminar at FIAN (Moscow). 
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